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Abstract 



The derivation of the Backlund transformations (BTs) is a standard problem of the theory 
of the integrable systems. Here, I discuss the equations describing the BTs for the Ablowitz- 
Ladik hierarchy (ALH), which have been already obtained by several authors. The main aim 
of this work is to solve these equations. This can be done in the framework of the so-called 
functional representation of the ALH, when an infinite number of the evolutionary equations 
are replaced, using the Miwa's shifts, with a few equations linking tau-functions with different 
arguments. It is shown that starting from these equations it is possible to obtain explicit 
solutions of the BT equations. In other words, the main result of this work is a presentation of 
the discrete BTs as a superposition of an infinite number of evolutionary flows of the hierarchy. 
These results are used to derive the superposition formulae for the BTs as well as pure soliton 
solutions. 



1 Introduction. 



The present paper is devoted to the Backlund transformations (BTs) for the Ablowitz-Ladik hier- 
archy (ALH) Since the discovery of the inverse scattering transform (1ST) these transforma- 
tions have become one of the most powerful tools for generating solutions for nonlinear integrable 
equations (one can find the discussion of various aspects of this subject in the monograph 3 ). 
The scheme of the solution generation can be briefly described as follows. One of the main points 
of the 1ST is the so-called zero-curvature representation (ZCR), when a nonlinear equation is pre- 
sented as a compatibility condition for some overdetermined auxiliary linear system. To this linear 
system one can apply the Darboux transform. This leads to some new solutions of the auxiliary 
problem and, consequently, of the initial nonlinear equation, which are determined (implicitly in 
general case) by the solutions we started with. In other words, the Darboux transformations for 
the auxiliary linear equations give us the algorithm to construct the Backlund transformations for 
the nonlinear one. This scheme was applied in 1970s to almost all known at that time integrable 
nonlinear equations (one can find the main results of these studies in the book |J and references 
therein) . This construction has been implemented by several authors for the ALH [SJ El El EI w h° 
obtained the systems describing the corresponding BTs. 

The main aim of this work is not to derive these equations but to try to solve them. To do this 
I will combine the standard algorithm based on the ZCR with another approach to the BTs, which 
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is associated mostly with the KP or 2D Toda lattice equations and which is based on rewriting 
the equation in question in terms of the so-called Miwa variables (see below). One can find a 
rather detailed description of this technique in the paper by Adler and van Moerbeke During 
a long time the 1ST and the Miwa's representation were used, in some sense, independently: some 
results have been derived using one approach, some using the other. One of the motives to write 
this work is to demonstrate that the combination of these two methods is a very useful strategy 
which gives possibility of solving some problems which are difficult to solve using only one of them. 
The BTs for the ALH discussed in this paper is a bright example of such problem. It turns out 
that combining both these approaches is rather fruitful and gives possibility of obtaining formally 
explicit expressions for the BTs which is the main result of this paper. 

The plan of the paper is as follows. After deriving the equations describing the BTs (section 
0) and presenting the necessary facts related to the functional representation of the ALH (section 
13) I will obtain the formal solution of the Backhand equations, i.e. present the explicit realization 
of the BTs for the ALH (section . These will be used to calculate the superposition of several 
BTs (sectional and then to rederive the AT-soliton solutions (section 0. 



2 Derivation of the Backlund equations. 

All equations of the ALH, which is an infinite set of ordinal differential-difference equations, can 
be presented as the compatibility condition for the linear system 

tfn+l = U n <f n (2.1) 

d^ n = V^* n (2.2) 

where dj stands for d/dzj, ZjS are an infinite set of variables (times), ^> n is a 2-column (or 2x2 
matrix), U n is given by 

U n = U n (X)=^ n (2.3) 

and Vn s are some 2x2 matrices whose elements are some polynomials in A, A -1 and depend 
on q n , r n , Q n ±i, r n ±x, ... In what follows we will not use their explicit form (one can find how 
to construct these matrices in the pioneering work pQ or, e.g., in the book jlO|L To provide the 
self-consistency of the system (|2.1|) . <|2.2|l the matrices V n have to satisfy the following equations: 

d j U n = V<$ 1 U n -U n V& (2.4) 

which, when rewritten in terms of q n s and r n s, constitute the ALH. 

From the viewpoint of the 1ST, the BTs are the transformations of the form 

= M n $„ (2.5) 

which are compatible with the n —> n + 1 shift. This means that \t„ should satisfy the equation 
similar to (|2.1|l with the matrix U n being replaced with some other matrix U n of the same structure, 
which leads to the following equation for the matrix M n : 

M n+1 U n - U n M n = (2.6) 

where 

Rewriting this matrix equation as a scalar system for the elements of M n , 
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M n = °" J 1 (2.8) 

one comes to the equations 

Aa„+i + q n b n+1 = \a n + r n c n (2.9) 

\~ 1 d n+ i+r n c n+1 = X^dn + qnbn (2-10) 

A _1 6„+i + r n a n+ i = \b n + r n d n (2-11) 

Ac„+i + q n d n +i = \~ 1 c n + q n a n . (2-12) 

A simple analysis of (|2.9|) ~ (|2.12|) leads to the conclusion that there exist two types of solutions of 
this system (compare, e.g., with the transformations of the first and second kinds in the paper 0). 

2.1 B -transformations. 

This type of transformations is determined by the following choice of the dependence of M n on A: 

M n =( X2 "± a ' an X Z n )- (2.13) 



\aq n --i 



System (|2.9j) (I2.12j) now becomes 



a' (a„+i - On) + a (q n r n+1 - q n -ir n ) = (2-14) 
ar n+ i + a' r n a n+ i - dr n = (2-15) 
tt?n-i + a'q n a n - dq n = 0. (2-16) 



By introducing the tau-functions of the ALH, 



q n = — , r„ = — i t„ = t„_it„ + i + p„cr„ (2.17) 



one can find a n , 

^~n— 1 



'n^~n— 1 

while the remaining equations can be rewritten as 



(2.18) 



ap n T n -i + a'p n -iT n - dr n p n -i = (2-19) 
aT„CT„_i + a'r n _iCT„ - da n T n -i = 0. (2.20) 

These equations are enough to meet (|2.6ll . However, the condition (t„) 2 = T n ^iT n+ i+p n a n imposes 
some additional restrictions on the transformed (tilded) tau-functions, which can be presented in 
the form 

a'pn<rn - d'r n T n + dT n+ iT n -i = (2-21) 

The system (I2.19|) - (|2.21|l completely determines this kind of transformations, which in what follows 
I will denote by the symbol IB instead of tilde, B / = /. 



3 



2.2 B -transformations. 

This type of transformations corresponds to the following choice of the M-matrix: 

In what follows all quantities related to this kind of transformations, which I will denote by 
the symbol B , will be marked with overbars, which does not mean complex conjugation. After 
calculating a n , a n = T n -\T n /r n T n -x one can derive equations similar to (|2.19[1 - (|2.21|) . 

d~T n p n _ 1 + d! Vn-iPn - ap n T n -i = (2.23) 
dcx n+ rf n + d! a n T n+ i - a r n+ ia n = (2.24) 
aT n+ i?„_i - a' T n T n + d! a n p n = 0. (2.25) 

Thus we have obtained equations describing two types of elementary Backlund transformations: 

{ap„(Br„_ 1 ) +a' p n -i (Bt„) - dr n (B p n _{) = 
aT n (Ma n - 1 )+a'T n - l (Ma n )-da n {MT n - 1 )=0 (2.26) 
a 1 Pn ( B a n ) - d! r„ ( B r„) + dr n+1 ( B r„_i) = 

and _ _ _ _ _ 

!d_T n+ i (M_p n ) + dl_ T n (Mpn+i) - ap n+ i (Bt„) =0 
da n+1 (Bt„) + J' (T„ (B -r n+1 ) - ar n+1 |l a„) = (2.27) 
a r n+1 ( B 7»_i) - a' r„ ( B r n ) + d! a n ( B p„) = 0. 

Writing equations l|2.26|l and (|2.27|l is a crucial step in constructing the BTs for the equations 
in question. Now one can, in principle, forget about the ZCR for the ALH which was our starting 
point and consider the bilinear systems l|2.26|l . (|2.27|l as the definitions of our transformations 
(which could be derived by other methods, say, the Hirota's bilinear approach). 

Usually from the practical viewpoint the BTs are used as follows. One takes simple solutions 
of the integrable system in question (the ALH in our case) and then tries to solve the Backlund 
equations (equations H2.26(l or (j2.27jl in our case) to obtain some more complicated ones (probably 
the first who applied BTs for generation of solutions of soliton equations were Seeger et al [TT]). 
It is possible to solve the Backlund equations explicitly if we start with some trivial solutions 
(e.g., constant or pure soliton), but to do it in general case is a task not less difficult than to 
solve the initial one. Nevertheless, it turns out (and it is the main point of this work) that it 
is possible to derive some (maybe formal) solutions for the Backlund system (equations (|2.26(l or 
(I2.27(l h This is difficult to do if one restricts oneself to some particular integrable equation (say, 
to the discrete nonlinear Schrodinger equation or discrete modified KdV equation, both of which 
belong to the ALH). At the same time, information contained in all equations of the hierarchy is 
sufficient to tackle this problem. We can construct a discrete (Backlund) flow of an infinite number 
of continuous ones. In other words, if we know that our tau-functions are solutions of the ALH (i.e. 
of the infinite set of compatible differential equations) , then we can use them to derive solutions of 
(|2.26|l and l|2.27(l . This will be done in section 01 after presenting some basic facts about the ALH 
in the next section. 

3 Ablowitz-Ladik hierarchy. 

Four our purposes we need the so-called functional representation of the ALH, which has been 
elaborated in |12l ITS*] , when the infinite set of differential-difference equations of the ALH is pre- 
sented as a few functional equations relating tau-functions of different (shifted) arguments. This 
way of rewriting a hierarchy using the Miwa's shifts (in the form of the so-called Fay's identities) 
is very convenient when one studies the hierarchy 'as a whole' especially in the case like ours when 
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we want to deal with all evolutionary flows simultaneously. The reader can find the derivation of 
this form of the ALH together with some formulae which will be used below in ^3]. However, to 
simplify the following calculations it seems fruitful to rewrite the results of |13| in a more general 
way. To this end it should be noted that some formulae of the previous section, as well as of > 
are related by the shifts p„ — > r„ — * a n . That is why I rewrite the tau- functions as 



Pn T n T n ~ T n 

and introduce a double infinite set of tau-functions t£ 
the relation written in l[2.17|) . 



in 



-- rl (3.1) 
±2, ±3, ... using a generalization of 



(C) 2 



m m i m — 1 m+1 
'n— 1 n+1 ' n n 



Thus I define 



and 



_m-2 



1 



.m+2 



(c +1 r 



m—l m — l 
~n-l T n+1 



m+1 m+1 
T n-1 T n+1 



for 



for 



2,3, 



-2,-3, 



(3.2) 
(3.3) 
(3.4) 



It can be shown that all equations of the ALH are compatible with l|3.3|) and (|3.4() and all equations 
of the extended version of the ALH can be obtained from the following ones: 

- 1 (E f O-C(E c C +1 ) (3.5) 
(E C C +1 )-C(E C C) (3.6) 



= 
= 



C-i(E<C+i) + 



i~ T n 



(the positive subhierarchy) and 



m-\-l\ 



t ~ m / IT? ~ m -r- 
CT„+1 (.H T n-l 



'n+1 



'E, 



jm — 1 



m+1 



m+n 



m+l\ 



(3.7) 
(3.8) 

(3.9) 
(3.10) 

(3.11) 

The extended version of the ALH l|3.5|) (|3.8|) was written in the paper ^1] by Sadakane who showed, 
using the free fermions approach, that the ALH arises from a reduction of the two-component 
Toda lattice hierarchy, discussed the ALH in the context of the affine Lie algebra that acts on the 
universal Grassmann manifold and clarified the universality of the ALH. 

In what follows, we also need some formulae from describing superposition of Miwa's shifts. 
They are collected in their extended form (i.e. in terms of the tau-functions of the extended ALH 
hierarchy, r™) in the appendix. 

Now we have all necessary to 'solve' equations (|2.26|l and l|2.27|l which can be rewritten as 



(the negative one). Here the symbols E^ and E„ stand for the Miwa's shifts 
E c f(z, z) = f(z + i[Q,z), E c f(z, z) = f(z, z + i[C]) 
which are defined for functions of an infinite number of variables, 

f(z,z) = f (Z1,Z 2 ,Z 3 , ■ ..Z 1 ,Z 2 ,Z 3 , ■■■), 

by 

f{z + a[£],z + 0[rj\) = f{z x + o& z 2 + a?/2, z 3 + a£ 3 /3, . . 

z x + fir], z 2 + (3r] 2 /2, z 3 + . . 



a C +1 ( 1 C+ 1 ) + a' r™ ( B C++ 1 ) - dr^ 1 ( B r™) = 



r m+n 

' n ) 
-m — l ( to> _m+l^ 



(BO+dC+i 



= 



and 



dr™^ 1 (Br™) + J'r™ +1 (] 



'n+1 



'n+1 



'n+1) 



'n+1 



r m+n 



a'C(Br™)+J'r™+ 1 



(3.12) 



(3.13) 
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4 Solution of the Backlund equations. 



As it was mentioned at the end of section |2 one can obtain, starting from the fact that r™ 1 solve 
the ALH equations, some particular solutions of the Backlund equations l|3.12[l and (|3.13|) . Indeed, 
by multiplying (|3.5|l by (d/a) m (d/d') n one can come to the conclusion that the quantity 

solves the first of equations 1)3. 12[) provided 

In a similar way, one can check that (|4.1J) also satisfies the second one. Thus, we have derived some 
BT which can be expressed in terms of the Miwa's shifts (i.e. the evolutionary flows) accompanied 
by the shifts of the site index n. Of course, the transform B^ is almost trivial and is not interesting 
from the viewpoint of applications. However, this transformation is not the only one and it is 
possible to derive three other similar transformations, involving both positive and negative shifts: 



d\ m fd 



d' 



E C C +1 (4.3) 

» (2) C = (5) (-^"Ec'C" 1 (4-4) 

Ei/fC (4.5) 



a 



d 

-) E^r^ 1 (4.6) 



with 



<~3 <«» 



which can be used to construct more rich ones, 

4 



.t^^/'iWt;, (4.8) 



fc=l 



It is an important moment. Each of B^r™ solves the bilinear ALH equations (|3.5[l 13.8fl . However 
their linear combination does not satisfy them automatically. Nevertheless, the integrable bilinear 
equations possess some peculiar features and it turns out that by imposing some restrictions on 
the functions in l|4.8|l one can make B t™ to solve the ALH equations. To do this, one needs 
to study the properties of the products of the Miwa's shifts and the trivial BTs B^^. 

From the superposition formulae for the Miwa's shifts presented in the appendix one can deduce 
that all B( fc ) satisfy the similar identities: 

aW(OC (e c b«c) = (EeO (b«c) +£ £ (Eer^O (bWc_ x ) (4.9) 

where 

A«(£) = a(C,0/C (4-10) 

A< 2 )(£) = c&C) (4.11) 

A (3) (C) = e(e,l/C) (4.12) 

A (4) (0 = fl(e,l/0- (4-13) 
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Thus, if one takes the functions uW depending on the variables zj according to the rule 

Eg _ AW(Q 

where A(£) is an arbitrary function satisfying the restriction A(0) = 1, then Br™ satisfies 



1 



A(0 E ? ( B O = — ( E £ O ( B O + e -7 ( E 5 r™ +1 ) ( B r™ : 



Now by straightforward algebra one can show that the transformation 



(4.14) 



(4.15) 



(4.16) 



is compatible with the evolutionary equations: if one substitutes in, e.g., (|3.5[) r™ with t™ (which 
is a solution of the Backlund equations (|3.12|) ') and Eg r™ with the r.h.s of H4.16[) then the resulting 
expression will be identically zero. 

Analogously, B t™ solve equations of the negative subhierarchy, if the dependence of the func- 
tions Uk on the variables Zj is given by 



where 



/,«(„) 
M (3) W 
M (4) W 



e(C^) 

s(C,»?) 

Ca(i/C,»7) 

c(i/C,v) 



and jti(??) is an arbitrary function with /i(0) = 1. 

Equations (|4.14|) and (|4.17|l can be easily solved by the ansatz 



(4.17) 

(4.18) 
(4.19) 
(4.20) 
(4.21) 



,(fc) 



(z, z) = exp < ^2 Ukj'< 



UJkjZj 



where 



and 



-i In 



(4.22) 



(4.23) 



(4.24) 



i.e., cjfej and Wjy are the Taylor coefficients for the logarithmic derivatives of the functions \^/\ 
and //W jii: 



(4.25) 
(4.26) 
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These formulae are nothing but the dispersion laws for the given boundary conditions (which 
determine the quantities a(£, r]), e(£, rj), etc, and hence all the functions and n^ k \r])). 

Similar algorithm can be developed for the B transformations. Starting from almost trivial 



r( 2 ) 



™(3) 



d 

a 
d 



y \ m 



-l/C '"+! 



with 



c 



one can construct the more general transformations 



k=l 



(4.27) 
(4.28) 
(4.29) 
(4.30) 

(4.31) 
(4.32) 



where the functions should be determined from the equations 

E^W _ A<*>(£) E^W _ /i^) (77) 



HO 



ith 



A (1) (0 
A( 2 )(0 



ff(e,c") 

e(£,C) 

c(C,l/C) 

Ca(VC,0 



M (3) W 
M (4) (l) 



c(»7, C) 

a(C,»7)/C 



(4.33) 



(4.34) 



and arbitrary A(£) and p,(rj) satisfying A(0) = /2(0) = 1. 

So, we have derived the main result of this paper: we have explicitly constructed the BTs for 
the ALH. This was done by combining the elementary transformations n — > n± 1, m — > rail, the 
Miwa's shifts and the multiplication by linear in n and m exponents. Each of these transformations 
is trivial and does not change the structure of solutions. However their combination can produce 
more rich BTs which give us possibility of obtaining more complicated solutions from simple 
ones. I will return to this point in the conclusion and discuss it using as an example the soliton- 
adding transformations of section |5| An important moment is that to obtain physically interesting 

transformations one has to use superposition of B( fc ) and B^ with different ks. Moreover, one 
cannot restrict oneself with k — 1,2 or k = 3, 4 only. Because of this it is impossible to obtain 
their infinitesimal versions by taking, say, the £ — > and £ — > limits. That is an illustration of 
the discrete character of B and B , contrary to the evolutionary operators E and E . 



5 Superposition of BTs. 

In the previous section we have derived the BTs, the transformations which commute with the 
n — > n+ 1 shifts as well as with the evolutionary flows. Now our aim is to ensure the commutativity 
of these transformations. It turns out that this problem is, in some sense, more difficult than the 
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previous ones. The commutativity with the n — > n + 1 shifts is due to the construction (|2.6() : 
all transformations H4.3j) - H4.6|) and l|4.27|) - [|4..WJI automatically satisfy (|2.26|) and (|2.27jl which was 
derived from H2.6JI ■ The commutativity with the evolutionary flows was achieved by proper choice 
of the dependence of the functions and vS k ^ (k = 1, 4) on Zj and Zj (j = 1, 2, ...) which is 
described by equations l|4.14l) . (|4.17l) and 1|4.33[) . 

Before proceeding further I have to introduce the following notation. Since now we will deal 

(k) 

with different BTs, which will be distinguished by the additional subscript. The symbol means 
the fcth elementary transform (one of (|4.3(l - (|4.6|l 'l constructed using the set of parameters a,j, a'j, 
dj and d'j. Consequently all the quantities £, the functions X^ k \ vP*', A, and others also will posses 
this index: Q — —a'jdj/ajd'j, X"- (£) = A(£; aj, a'j, dj, dj) etc. Finally Mj is the BT corresponding 

(k) (k) 

to the superposition of with coefficients Uj . The similar notation will be used for the BTs 
of the B-type. 

Now we can return to the superposition of BTs and I would like to start with an important 
remark. It should be noted that the commutativity of evolutionary flows and BTs has been done 
for B and B being taken as combinations of four elementary transformations, see ()4.8|l and (|4.32|) . 
However, as will be shown below, if we want B and B to commute between themselves, we have 
to restrict ourselves to combinations of only two elementary BTs. Indeed, if we apply the to 
the B 2 -transformed tau- function B 2 r™ , we get 

••?(*;) r m J (M[ k) M 2 T^ = — (v u "r"< \ i v ,-»<- i ^ 'h_ 

where 

(5.2) 

(5.3) 
(5.4) 
(5.5) 

A simple analysis leads to the conclusion that it impossible to chose the functions A 2 (£) and ^{ff) 
satisfying A 2 (0) = /i 2 (0) = 1 in a way to met the condition 

= 0$ = P { i2 = P ( i2 (5-6) 

and to ensure the skew-symmetry of the left-hand side of (|5.1|l with respect to the interchange of 
the subscripts 1 and 2. However there is no problem, if one restricts oneself to only two elementary 
transformations. For example, if 

Mj =uf ) mf ) +uf ) mf (5.7) 

(this is the situation which corresponds to the bright soliton case, see also section EJ) then the 
choice Aj(£) = 1, fJtj(rj) = 1 + 1] dj/aj together with the restriction a'j = dj leads to 

Pi2 C _1 (Bi2 = -(BiT™) (B2T™- 1 ) - * (Brr™- 1 ) (B 2 r™) (5.8) 
a 2 a\ 

with (3\2 = d 2 /a 2 — d\ja\. Hereafter, the multi index of B indicates superposition of a few BTs 

ly...C= BiB,- ...r™ (5.9) 







^ = 


^A 2 (d) 

0-2 


0[f = 


d 2 1 - C1/C2 
a 2 A 2 (Ci) 


= 


Ci Ml/Ci) 


f® = 


4 1-C2/C1 

a 2 M2(l/Ci) 



9 



(after we have ensured the commutativity of different BTs, this notation has sense, because the 
result does not depend on the order of the application of BTs) . 
In a similar way, for 



,(3) 1»( 3 ) 



(5.10) 

■rfa'j/aj 



(this is what one has in the dark soliton case), the choice = 1 + £ctj/a'j, = 

together with the restriction dj — dj leads again to l|5.8|l . with (3± 2 — o! 2 ja 2 — a[/ai. 

To conclude this analysis I would like to note that we are at the point where the boundary 
conditions play a crucial role. They entered earlier specifying the dispersion laws through the 
quantities a(£, 77), etc (see the previous section). But now their effect is a 'qualitative' one. They 
select the structure of the BTs and impose some restrictions on their parameters. 

In the following consideration I will not write explicitly the values of the coefficients like f3jk 
(though we already know that in the most relevant cases /3jk = Cj — Cfc) because the 'physical' 
quantities, such as q n , r n of the ALH in the original setting, depend on ratios of the tau-functions. 
So, the constant (both with respect to the indices n, m and the times) multipliers ((3jk and similar 
ones which appear below) are not crucial. 

By applying the Backlund equations l|3.12fl and i|3.13[l to the superposition formula (|5.8|) derived 
above, one can obtain a large number of the equivalent ones. Some of them are given by 



/?12C +1 (Bl2C) 

/3i2t:- 1 (Bi 2 C +1 ) 



C) 



^0 

— — (Bi T% 

ai a 2 



;C +1 )-^(Bir^ 1 )( B 2C) 
ai 



n+l) 



d[d 2 

— — i^lT, 

Oi a 2 



(b 1 <; 1 )(b 2 c). 



(5.11) 
(5.12) 



In a similar way, one can derive formulae describing the product of the B transformations 



^i2C +1 (B 12 rr) 
&C +1 (»i2r;) 



IT, 



l )0 



«2 

d 2 
d 2 K 



fll fl 2 (m J, 

di d 2 



m+1) 



«1T, 



m+1 



l )0 



n+l) 



ai 
" di 

^(iiC +1 )(i a r;- 



)(B 2 r™ 1 )-H(BiC +1 )(B 2 < 



(5.13) 
(5.14) 
(5.15) 



and the 



products, 



712 C(I 



712^(8! B 2 r™ 



^C+^Ba B 2 r™ 



(Bi< 

d[d 2( 
d\ d 2 

— — (Bit, 
di a 2 



d\ a 2 



m+l\ 



»2T„ 



m+l\ I 



n t„. 



4 , 
a 2 



s 1 o(B 2 c + - 1 1 ) 



n+l 



«2 T„ 



(5.16) 
(5.17) 
(5.18) 



Here, again, j3\ 2 and 712 are some constants which depend on the boundary conditions and their 
explicit form should be established after specifying the class of solutions we are dealing with. 

After we have obtained explicit formulae for the superposition of two BTs it is easy to generalize 
them to an arbitrary number of BTs involved by noting their determinant structure. For example, 
equation l|5.11|l can be rewritten as 



6 12 T n 



1 

/3l2 T, 



1 



det 



n+l 



ai V 



;b 2 o 

(Bar™.. 



(5.19) 



and generalized as 



l [N] 



T[N] 



T-riV-l 

rL- t. 



det 



n+j 



(B fc r 5 



n+j -I) 



(5.20) 
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where the traditional designation T™[iV] is used for the N times transformed tau-function, 

r™[N] = M 1 ...M N r™ (5.21) 

and the constant T[N] is given by 



T[N] 



n 

l<j<k<N 



(5.22) 



A proof of this formula can be given as follows. Subtracting from the jth column, j = 1, N — 1 
the (j + l)th one multiplied by the factor (ajsi/a' N ) ( Bjy T ,7+j-i/ ®JV T n+j) one can make all elements 
of the Nth row equal to zero, except the element at the (N, N)th place. In such a way the Nth 
order determinant at the right-hand side of (|5.2U|I is reduced to the one of order N — 1. The 



elements of the resulting determinant can be shown to be proportional to B jvfc t™ 



So, the 



larger determinant can be presented as the result of the application of B^r to the smaller one. 
Repeating this step one comes to (I5.20f) . 

In a similar way, one can obtain other formulae describing the superposition of B -transfor- 
mations, B -transformations as well as the 'mixed' ones for the superposition of both B - and 
B -transformations. 

The superposition formula for the join action of both B - and B -transformations can be written 

as 

T[N, N] 



11 7 =1 'n4 



n+J 



detMi 



Here 



t™[N,W] = Bi ... Bjy Bi...%r™, 
the numbers M and M, M + M = N + N, are given by 

M = N -m + m, M = N + m - fi 

and M% is the (N + ~N) x (N + N) matrix 

A»[M,N\ B£[M,N] 



Mi 



C%{M,N] V%{M,N] 



(5.23) 
(5.24) 
(5.25) 

(5.26) 



(the numbers in the square brackets determine the sizes of the matrices that form the matrix 
A£[M, N] is a M x N matrix, etc). The elements of the matrices A, B, C and T> are given by 



K[M,N] hK 

B£[M,7?\ 3 ,* 
C%[M,N]j, K 



a K 

Ok 

d K 



8 T P 



.4 



d f , 

3-1 



J-l 



u+l 
T n+j-l 



In all above formulae 

.7=1, 



J=1,...,M K=1,...,N K = l,. 



(5.27) 
(5.28) 
(5.29) 
(5.30) 

(5.31) 
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and the factor T[N, N] is given by 



T[N, N] 



(-) 



(m-fi){N+N) 



/3[N] j[N, N] P[N] 




m — fi 



with 



0[n] = n 

1<3<K<N 



N N 

j[N,N]=l[l[ 7 y l 

3=1 K=l 



0[N] 



n j*i 

l<j<ii<N 



(5.32) 



(5.33) 



To illustrate this general result let us consider the simplest case when is a 2 x 2 matrix. All 
binary superposition formulae presented in this section can be obtained from l|5.23|l . The following 
table contains the values of the parameters N, N and p which lead to l|5.8|l . (|5.11|) - (|5.18|l . 



N 


N 


fj, 


equation 






m — 1 




2 





771 


(|5~TTf 






m- 2 


ifB~T2) 






777, + 1 







2 


m 


ii5~n 






777 +2 


(|5~TB] 






777 


||5jJ| 


1 


1 


777 —1 


(|o~T7ji 






777 +1 





6 Solitons of the ALH. 

In this section I will discuss two types of soliton solutions: the bright solitons (which correspond to 
the case of the zero boundary conditions) and the dark ones (the case of the so-called finite-density 
boundary conditions). They will be obtained by applying the A^-fold (2iV-fold) BTs to the vacuum 
solutions. It should be noted that the formulae of the previous section immediately give us the 
structure of iV-soliton solutions. So, we have only to write the dispersion laws and to solve the 
question of commutativity of different BTs. 

Contrary to traditional IST-based approach I will distinguish different types of solitons not by 
writing the boundary conditions for the functions q n and r„, but starting from the corresponding 
vacuum solutions (0-soliton ones) directly in terms of the tau- functions r™. 



6.1 Bright solitons of the ALH. 

The vacuum tau-function corresponding to the bright soliton case is given by 

_m s 

T n — °rn.O 

or, if to return to the original tau- functions, r n , a n and p n , by 



Tn — 1 , On — Pn 

Applying the elementary BTs l|4.3|l - H4.6|l . one gets 



0. 



BW T n = (d/d') n 

B( 2 ) T n = 

' T„ = (d/d'Y 

T„ = 



(1) 



= 



{2) o n = {-a/a'T 
(3) a„ = 



o n = {-a/ a') 



An 



Pn 
Pn 
Pn 
Pn 



(6.1) 



(6.2) 



(6.3) 
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It is easy to note that 

lW=B( 3 ), B^=B( 4 ) (6.4) 
which means that general BT, B , can be taken as a combination of two elementary one, say, as 

M =u^M^ +u^M ( - 3 '> (6.5) 

and to ensure the commutativity it is sufficient to impose the restriction 

a' = d! (6.6) 

and to choose the functions A(£) and n{rj) as follows: 

A(0 = 1, fx(v) = l-vC, C = -- (6.7) 

a 

which leads to equations l|5.8|) . (|5 . 1 1|> and (|5.12() with 

/3i 2 =Ci-C 2 . (6.8) 
Calculating the coefficients that appear in the superposition of the Miwa's shifts 

a(Z,v) = a(£,ri) = £ - V, e(t,rj) = 1- (6.9) 

and 

c(Z,v)=£(Z,r))=9(Z,v)=l (6-10) 

which is easy to do by the substitution of H6.1|l in the formulae presented in the appendix, one 

comes to the following equations describing the dependence of the functions on the variables 
of the ALH: 



Eg u (Z> Eg u {Z > i 

Eg . cjr E, u (3) 1 

Similar analysis can be applied to the B -transformations. After calculating the elementary 
transformations 

it /ii /u 



(6.12) 



: (-d/d') n 

(6.13) 
(6.14) 

which leads to the restriction 

a' = <? (6.15) 



B (1 V„ 


= (a/aT 


B^(7 n = 





B (1) p« 


B (2) r„ 


= 











B cr„ = 


B (2) Pu 


B (3) r„ 


= (a/aT 











B er„ = 


B (3 Vn 


B (4) r„ 


= 











B <7 rl = 


B (4) p„ 


and noting that 














B (1) =B (3) , ' 


B (2) = I 




one can choose 




I = u 2 B (2) H 


- m 3 B (3) 





















and 



A(0=1, = 1 - »?/C, C = -J (6.16) 
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Further, one can get (3 12 — Ci ~ C2 and the following equations: 



= 1 
= 1 



1 



(6.17) 



An important note: both B - and B -transformations, taken alone, destroy the symmetry be- 
tween q n and r n . Say, B q n — const - C _n j while lr„ = 0, which means that it is impossible to 
make the transformed functions to satisfy the 'physical' involution: q n = — r* (the star stands 
for the complex conjugation). The same is valid for the B -transformations. Hence, to obtain 
'physical' solutions in the bright soliton case we have to use the B - and B -transformations only 
in combination, which gives possibility of preserving the involution. Say, to obtain the iV-soliton 
solution we have to use the superposition of ./V B - and N B -transformations (with imposing some 
relations on parameters of the former and the later). Setting 



N = N 



and noting that 



and 



T„ = — »„ T 



C„ = r J» K (J 



rr BftT 



one can rewrite expression (|5.23[) for the superposition of BTs, in the /1 = case, as 



1 [N, N] = (-)^ m(m+1) C,CC 
x det 



n-3 , 



i« C ) 



TO 



3=1, 
3=1, 



,N—m 
,N+m 
,...,N 



(6.18) 
(6.19) 
(6.20) 

(6.21) 



where the quantities £*, and C* are given by 



= ( _ ) 1 r \ 



I\k=1 C« 

ris=i Ck 




and 



(-) 



|JV(JV-1) 



/3[7V] 7 [7V, N]/3[iV] 
That leads to the following formula for the TV-soliton solution: 

A (i) 



g„[JV,iV] = -C» 



A (0) ' 



where 



Ai B) =Ai ro) (z,z) = dot 



r n [N,N]=Q 



Rr (z, z) C s 



ci- 1 



(-1) 



A (o) 



ra-j 



j—l,...,N—m 
3=1,.. -,N+m 

K,R=1,...,N 



(6.22) 



(6.23) 



(6.24) 



(6.25) 
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The quantities Q K (z, z) and R R (z, z), which describe the dependence on the ALH-coordinates z 



and Zj, 



(2) -(2) 
}k{z,z) = -£rr, i? s (z,z) = -^r, (6.26) 



can be easily calculated from l|6.11|l and (|6.17|) using the identity 



ln(l -*) = -£- (6.27) 



and are given by 



Q K (z,z) = O:: < xj> I ' ' j (6.28) 

ife(M) = 4 0) ex P | i£(%Ci + %C^)J (6-29) 

where and i?^ are some arbitrary constants (these expressions are nothing but the well-known 
dispersion laws for the bright solitons). 

6.2 Dark solitons of the ALH. 

The vacuum tau-function corresponding to the so-called finite-density boundary conditions (the 
dark soliton case) is given by 

r™ = a^p^u m v n (6.30) 

where the constants a and (3 are related by a + (3 = 1, while u and v are some functions of 
the ALH-coordinates. Contrary to the bright solitons which are two parametrical, a dark soliton 
depends only on one parameter. From the viewpoint of the BTs this difference manifests itself in 
the following facts: (1) to construct A^-soliton solutions we need to apply N BTs to the vacuum 
solution (and not 2N, as was in the case bright solitons) and (2) the B and B transformations 
lead to similar results. So, in what follows the A^-dark soliton solutions will be obtained as a 
superposition of the B transformations using (|5.2U|) . Thus, we have to solve now the following 
three problems: we have to ensure the commutativity of our BTs, to calculate the dispersion laws 
and to satisfy the 'physical' involution condition r n = 5*. As in the previous section, the first 
problem can be solved rather quickly. The case is that of four elementary BTs H4.3H ~ H4.6H oru y 
two, say BW and B< 3 \ are independent while B< 2 ) = B^ 3 ) and B< 4 ) = B^. So, one can define the 
B transformation as 

B = U «B« W 3) B (3) . (6.31) 
Then it is easy to show that they will commute, if we set 

d=d' (6.32) 

and choose 

A(0 - 1 - C/C, l*(v) = 1 - <V- (6.33) 
This leads to the superposition formula (|5.8|l with /3i 2 = Ci — C2- 
The dependence of w 1 ' 3 ' on the ALH times is given by 



E ; u (1) 1 E,, u (1) 1 

Eg u (3) 1 E v u {3) 1 



(6.34) 
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where c, c and g are the coefficients which appear in formulae describing the superposition of the 
Miwa's shifts E and E (see appendix). Their explicit form, as well as the explicit form of the 
functions u and v from H6.30|) can obtained by substitution of the vacuum tau-function l|6.3()|) into 
the corresponding equation. This gives 



g(£,v) 



and 



E ( u 



Ec u 



1 — a^rj 



E ( u 



( E v u^ 


) +/3 


V u , 




E^ u 




u 



E v v 



E c v 



R(0 - C + 1 

2 

W) + C + 1 



where 



R 2 (() = l + 2(a-f3)( + C 2 



(6.35) 
(6.36) 

(6.37) 
(6.38) 



The last equations contain all necessary to obtain the dispersion laws that we need. However the 
presence of radicals -R(C) complicates the following analysis, so, I would like to rewrite them using 
the parameterization 

sin 29(C) 



c 



sin [2w - 20(C)] 



where the constant parameter lo is given by 



a — cos lo, 



In terms of 0(C), 



and 



R(0 



(3 = sin 2 lo. 
sin 2lo 



sin [2w - 20(C)] 
cos [0(C) - 6>(r?)] 



and we can present B as 



with 



cos 0(C) cos 6(rf) 
sin lo sin [lo — 0(C) — 9(v)] 
sin [w - 0(C)] sin [lo - 0(? ? )] 

sin w sin(cij — 0) 




cos 
sin w cos(w — 0) 



sin 



(6.39) 

(6.40) 
(6.41) 

(6.42) 
(6.43) 

(6.44) 

(6.45) 
(6.46) 



(the functions and u/ 3 ' are, up to some unimportant constants, to' 1 ' and u^). 

From these formulae, one can deduce that the condition q n = r* can be satisfied by choosing 
d/a = cos lo and introducing the real parameters 7fc, describing the fcth BT, by 0fe — lo/2 + ijk, or 



sin (lo + 2ij k ) 
sin (lo — 2ijk) 



(6.47) 
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which reproduces the already known fact that the eigenvalues of the scattering problem, corre- 
sponding to the dark solitons are located on the ark of the unit circle |Cfc| = 1, | arg£fc| < 7r — 2uj 
(in our designations). Now the BTs Bfc can be written as 



(6.48) 



Here, the constant b k is given by b k — y vu^u^ sin — ijk) / cos + ij k ) (I write this ex- 
pression only for the sake of completeness, because b k does not appear in the final formulae) while 
other parameters (more essential) are given by 

■& k — argsin (to - 2i^ k ) , Xk = hi sin u + i<& k (6.49) 

and 

h k = tan + ijk^j , (f> k = argtan (~ + ij k ) . (6.50) 

The dependence on the ALH-coordinates Zj and Zj is 'hidden' in the functions w k — w k (z,z) 

(which are related to the old ones, u k 1,3 \ by w k oc \J vu k ^ /u^). This dependence is given by 
equations (|6.34() which lead to 



W k = W { k °> exp 1 ^2 (P-kaZa + /4a z a) 



(6.51) 



where w, are arbitrary constants and fx k are complex conjugate to fi ka while the last ones are 



determined by 



JV^ C Q \ ^n [f+t 7 «- 0(C)] 



tan 



(6.52) 

, .f + »7kJ 

Applying the £^ to the logarithm of the last equation one can conclude that the quantities \i ka 
are the coefficients of the following Taylor series: 



% sin 2uj £ 



2 sin26> K R(() 1 - C exp (2iz? K ) 



(6.53) 



Now we have all necessary to write the final formula for the A-dark soliton solution of the ALH: 



q n [N] = q vac exp (i ^ <d K 



A (+) 

7^ 



r n [N)=q* n [N] 



with 



A (+) 



det 



det 



K^- 1 + u;-^-"-^ 1 } exp (-i§ K j) 

j,k=l,...,N 

{e^w^- 1 + e-^w- 1 h- n -=+ 1 } exp (—i$ K j) 



j,k=l,...,N 



Here g„ 



a m is given by 



9«ac(2, z) = Va exp < i A a (z a + z Q ) 



I. o=l 

where A a are the coefficients of the Taylor series 

fw-'S-t- 1 

(which follows from l|6.37[) ). 



2i?(C) 



(6.54) 

(6.55) 
(6.56) 

(6.57) 

(6.58) 
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7 Conclusion. 



The main idea of this work was to analyse the equations describing the BTs for the ALH, using 
the so-called functional representation, when an infinite number of the evolutionary equations are 
replaced using the Miwa's shifts with a few difference equations. These equations contain much 
more information than any of the equations of the hierarchy taken alone and this gave us the 
possibility of solving the BT equations. The main result of this work is explicit (though formal) 
solution of the BT equations, i.e. explicit form of the BTs. The fact that this can be done is not 
surprising and is already known from the general Sato's theory. For example, in [2] one can find 
rather comprehensive discussion of how this method works in the context of the KP equation. The 
solutions of the Backlund equations obtained in this paper were constructed of four elementary 
transformations: n — > n ± 1, m — > m ± 1, Miwa's shifts and multiplication by the factor of the 
form a m b n . Each of these transformations is almost trivial. For example, if we apply them to 
the dark-soliton vacuum l|tj.30fl . they coincide and do not change the quantities q n and r n . In the 
case of the bright solitons the situation is even worse. Say, the transformation m — > m ± 1 (which 
is one of the 'elementary BTs' of ^1]) to the vacuum tau- function l|tj.l|) he will obtain that new 
tau-functions are given by a n = 1 and r n = p n — which means that both q n and r n are undefined: 
q n — 1/0, r„ = 0/0. The key ingredient of the BTs is the linear superposition of different B( fc ) 

and B , the superposition which is natural for the approach based on the 1ST technique 0]. 
Taking as an example the results of sectional one could note that the 'soliton' exponents (|6.28|) 
and H6.29J) appear when we solve equations determining the coefficients and of the linear 
combinations (|6.5|) and (|6.14(l . 

Another important question is the following one: are the BTs of this paper the standard 
Darboux-Backlund transformations or something different? As it was said in the introduction, in 
the modern theory of the BTs, based on the 1ST approach, the BTs appear as the results of the 
Darboux transformations applied to the auxiliary linear problems. In this paper, we have used 
another approach (trying to minimize the use of the auxiliary - intermediate, in some sense - 
problem) and after having rederived the Backlund equations did not return to (|2.1|) . Nevertheless, 
I would like to state here, without presenting explicit formulae describing the action of our BTs 
at the linear level, that BT discussed in the present work are indeed the well-known Darboux- 
Backlund transformations of the soliton theory. As an illustration to this fact, one can consider 
the results of section HO where the BTs were acting as soliton-adding transformations, which is the 
most typical example of the Darboux-Backlund scheme. 

Finally, it should be noted that approach of the presented work (to use functional representation 
instead of the auxiliary linear problem) is not so sensitive to the boundary conditions than one 
based on the 1ST and can be more easily modified to the case of other classes of solutions of the 
ALH (not only soliton ones but, e.g., quasiperiodic, Wronskian, and other) without necessity to 
elaborate from scratch the corresponding inverse scattering scheme. 

Acknowledgements. 

I am grateful to the referees of this paper for careful reading the manuscript, correcting misprints, 
adding important references and making suggestions on improving the presentation of the results. 

This work has been supported by grants: BFM2003-02832 (Ministerio de Educacion y Ciencia, 
Spain), PAI-05-001 (Consejeria de Educacion y Ciencia de la Junta de Comunidades de Castilla-La 
Mancha, Spain). 

Appendix: superposition of Miwa's shifts 

Here, I would like to present formulae describing superposition of Miwa's shifts. Some of them 
were derived in other can be obtained from those of [13] by simple algebra. 
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The superposition of the positive Miwa's shifts E can be described by 

a(Z,n) C(E«E„C+i) - e(E S C +1 )(E,r:)-t ) (%C)(E,r n n ' +1 ) (A.l) 

a(£, 77) r™ (EfE, = £ (E ? C+ 1 ) (E, C) - 7? (E 4 C) (E, C +1 ) (A.2) 

a(£, 77) r™ (EeE, r^ 1 ) = £ (E £ C-+ 1 ) (E, C) - 77 (E 4 C) (E„ C++ 1 ) (A.3) 

where a(£, 77) is a constant (with respect to the indices m and n) which should be determined form 
the boundary conditions. Note that only two of these equations are independent: the third one 
can be obtained by applying (|3.2|) . In the same manner one can produce a large number of similar 
formulae. Another type of superposition formulae for the positive Miwa's shifts can be presented 
as follows: 

c($, 7?) (Ee C) (E, C) = C-! (E«E, C+i) + C" 1 (E S E, C +1 ) (A.4) 
c(£, 77) (E S C) (E„ C+i) = r n ™ (EeE, + £ C" 1 (EcE, C-+ 1 ) (A.5) 

c(£, 77) (E 4 (E, C) = C (E?E, - ^7 C-i (E ? E, C++ 1 ) (A.6) 

where c(£, 77) = (f - ?/)/a(f , 77). _ 

The superposition of the negative Miwa's shifts E can be described by 

afo 77) c +1 C) - e (E c C) (E„ C+i) - 77 (E 6 r™ +1 ) (E„ C) (A.7) 
a(£, 77) c (I € I„ c +1 ) = e (e 5 C +1 ) (1, c) - 77 (Ej C) (I, C +1 ) (A.8) 

5(C, 77) C +1 (I € I„ C+ 1 ) = £ (E ? r™ +1 ) (1, r™ x ) - 77 (E> r™ (E, C +1 ) (A.9) 

or by 

c(£ 77) (E ? C) (I„ C) = r™ +1 (Iff, C-i) + C" 1 (E e E r; (A.10) 

c(e, 77) (E 4 C) (E, C+i) = C+i (I ? I r) C) + 77 C+1 1 (E C E ?; C +1 ) (A.ll) 

g(€, 77) (I ? c) (e„ C +1 ) = C (%!, c +1 ) - € C+i (e s I, C+1 1 ) (A.12) 

where, again, c(£, 77) = (£ - ??)M£, v)- _ 

Finally, for the mixed products of E and E shifts one can get formulae such as 



and 



e (e,r 7 )C(E€E 7) C) 


= (E ? C)(E, ; C)~ 


^(E £ C+i) (E,C-i) 


(A.13) 


e(e,77)C- 1 (E«E ?? C +1 ) 


= (E 5 C)(E,C)- 


(E ? C+i)(E,C-i) 


(A.14) 


g(Z, V ) (E^o^c) 


= C(E«E r) C)" 


£77 C" 1 (E?E,C +1 ) 


(A.15) 


g(t, V ) (E 5 C+i) (E„C-i) 


= C(E«E r) C)- 


C" 1 (E e E,C +1 ) 


(A.16) 


g&r,) (EeC+i^^C) 


= c+iMe^c) 


-^C^^C+i) 


(A.17) 


gfor,) (EeC+i^^C) 


= C(E«E J) <++ 1 ) 


-r/C+i (EeE,c +1 ) 


(A.18) 



with e(£, r))g(£, 77) = 1 - £77. 
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